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Results of the first principles study on the electronic structure and magnetism of the supercon- 
ducting weak ferromagnet Y4C03, are presented. Using the full potential Korringa-Kohn-Rostoker 
(FP-KKR) method, densities of states, dispersion curves and magnetic moments were calculated 
for quasi-ordered structural model of the compound in the framework of the local density approxi- 
mation. Spin-polarized KKR calculations confirm that weak ferromagnetic properties of Y4C03 can 
be attributed to only one cobalt atom located on (2b) site in the unit cell, while other twenty Co 
and Y atoms acts as a diamagnetic environment. Moreover, the magnetic Co atoms form a quasi- 
one-dimensional chains along 2 direction. The magnitude of Co(2b) magnetic moment (0.55 ^s) 
markedly overestimates the experimental value (0.23 /is), which suggests the importance of spin 
fluctuations in this system. Calculated distribution of spin magnetization in the unit cell provides 
a background for discussion of the coexistence of ferromagnetism and superconductivity in Y4C03. 
Finally, the effect of pressure on magnetism is discussed and compared with experimental data, also 
supporting weak ferromagnetic behaviors in the system. 

PACS numbers: 74.25. Jb,74.25.Ha,75.10.Lp 



I. INTRODUCTION 

Unusual properties of Y4C03, that is the coexistence of 
superconductivity and ferromagnetism, were observed 30 
years ago.^ The superconducting and ferromagnetic crit- 
ical temperatures are about Tg ^ 2.5 K and Tc 4.5 K, 
respectively. After this finding, the system was inten- 
sively studied using different techniques, since it was 
the first example where superconductivity coexisted with 
weak itinerant ferromagnetism, with both phenomena 
driven most likely by the d-band electrons. To our best 
knowledge, Y4C03 is up to now a unique compound 
containing transition-metal elements only, where coex- 
istence of ferromagnetism and superconductivity occurs, 
since other systems always contain /-like elements, e.g. 
UGe2 (Ref. 2), URhGe (Ref. ^), UCoGe (Ref. 4) or re- 
cently discovered P-doped EuFe2As2 (Ref. 5). More- 
over, Y4C03 was likely the first example, where super- 
conductivity occurred below a transition to ferromag- 
netism, with both phenomena coexisted at least in the 
range between 1 K and 2.5 This was another dif- 

ference to previously known magnetic superconductors 
e.g. ErRh4B4 or HoMogSg (see, Refs. Isl-flOl). where mag- 
netism appeared below T,, suppressing superconductiv- 
ity. 

Although Y4C03 has been known since 1980, the first 
principles study on its magnetism (e.g. calculations of 
magnetic moments) has not been presented until now. 
We have found two papers reporting electronic structure 
of this system from ah initio calculations. Interesting 
discussion of electronic structure of related compound 
Y9C07 based on non-spin-polarizcd computations can be 
found in Ref. illi (for connections between Y4C03 and 
YgCo7, see below). More recent paper, Ref. [ij, briefly 
reports electronic structure calculations for Y4C03, also 



in non-magnetic state. 

Noteworthy, one may notice a renewed interest on this 
system due to very recent experimental investigations^^ 

In this paper the results of band structure calculations 




FIG. 1: (Color online) (a) Unit cell of Y4C03 used in the 
calculations. The edges of the unit cell, i.e. (2b) sites, are 
occupied by cobalt or vacancy (in red), on average every sec- 
ond position is occupied by cobalt, called here Co(3). Atomic 
positions are presented in Table HI (b) The body of the unit 
cell is formed by the trigonal prisms of Y(l) (yellow), Y(2) 
(green), Co(l) (violet) and Co(2) (blue), the black lines and 
labels denotes atoms in the z = 0.75 plane, while white de- 
notes z — 0.25 plane, (c) The in-plane triangles of Y(l) atoms 
form channels along z direction, surrounding the Co(3) atoms 
chains. 
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for Y4C03 are presented. We attempt to answer some 
open questions, e.g. whether the magnetism is an intrin- 
sic property of the system or an effect of crystal imperfec- 
tions, -ii or what is the microscopic reason for magnetism 
occm'rence in Y4C03. The imphcations of our electronic 
structure calculations on the current model of the coexis- 
tence of ferromagnetism and superconductivity in Y4C03 
are also discussed. 



A. Review of the experimental data 

In this section we summarize shortly crystal structure 
and magnetism data, available for Y4C03. For a more 
complete review, see e.g. Refs. H0. The crystal struc- 
ture of the Y-Co system near 1:1.3 stoichiometry is rather 
complex, and was oryginally described in the hexago- 
nal unit cell of the H04C03 type^'* (Fig- [II- Unit cell 
can take 22 atoms, being distributed over three inequiv- 
alent Co sites and two inequivalent Y sites. All positions 
are gathered in Table HI As the unit cell includes three 
formula units, for the Y:Co stoichiometry equals to 4:3, 
the Co(2b) sites are half-filled (50%) and the number of 
atoms in the unit cell is equal to 21. Thus, in this crys- 
tallographic model, Y4C03 can not be regarded as an 
ordered compound, but as an disordered alloy with (2b) 
sites occupied randomly by cobalt (Co atom on (2b) site 
will be called Co(3)) and vacancies (Vac). Thus, the crys- 
tal cell of this alloy better corresponds to the formula 
Y12C08+2X with X = 0.50 than to the formula Y4C03, 
suggesting an ordered system. The unit cell (Fig. [1]) con- 
sists of the 'separated' Co/Vac chains forming the cell 
edges and the Y(l)-Y(2)-Co(l)-Co(2) trigonal prisms fill- 
ing most of the unit cell. The partial occupation of the 
Co(3) site is likely connected with a short distance be- 
tween the neighboring (2b) positions ('^2 A). 

In 80-ties, after revealing the coexistence of ferromag- 
netism and superconductivity in Y4C03, the H04C03- 
type unit cell appeared to be only an approximation for 
the real structure of the system. There was a debate 
whether this system should be called Y9C07 (instead 
of Y4C03), since it was found, ^'^^^^ that Y9C07 (richer 
in cobalt) exhibited better superconducting properties 
(^ 0.5 K higher Tg) than Y4C03. It was even suggested, 
that Y4C03 as a single phase compound might not ex- 
ist Employing the initial notation, the 9:7 stoichiome- 
try compound corresponds to Y4C03.11, but one should 
bear in mind that for Y9C07 an ordered model of struc- 
ture was suggested^^ with the unit cell tripled along the 
c-axis. In such a superstructure, the Co(3) atoms posi- 
tions are adapted in such way to make the Y-Co system 
perfectly ordered for 9:7 composition (100% filling of all 
sites, see e.g. Refs. [TlllTsI for details). The resulting crys- 
tal structure exhibits reduced symmetry, trigonal instead 
of hexagonalj^ However, the ordered model of the super- 
structure was not confirmed by X-ray measurements on 
single crystal samples.^"* Moreover, it was also reportedl"* 
that the Co (3) atoms did not occupy the fixed lattice po- 



TABLE I: Atomic positions of the Y4C03 compound (type 
H04C03, space group No. 176, P63/m). Lattice parameters 
areji^ a = 11.527 A, c = 4.052 A. 



Atom site 



coordinates 



Y(l) (6h) {x,y,\), {y,x-y,\) {y-x,x,\) 0.7543 0.9791 



Y(2) (6h) 

Co(l) (6h) 

Co(2) (2d) 

Co(3) (2b)" 



(2;,y, I) {y,y-x, |), {x-y,x, |) 

as above 
as above 

V3' 3' 4,'^ V31 31 4,) 

(0,0,0), (0,0, i) 



0.1360 0.5150 
0.4415 0.1578 



"In our calculations (2b) positions are split into (la) (0, 0, 0) and 
(lb) (0,0, ^), being occupied by Co and vacancy, respectively (see 
text). 



sitions, and were rather arranged in chains with undeter- 
mined period. Thus, the basic unit cell shown in Fig. [1] 
can be treated as the best established and simplest ap- 
proximation of the real structure both for Y4C03 and 
Y9C07, as far as the Co(3) atoms location is concerned. 
The Y12C08+22: model represents the Y4C03 and Y9C07 
compounds, when the Co(3) concentration reaches the 
values of a; = 0.50 and x — 0.66, respectively. However, 
the question on the real unit cell of Y4C03 and Y9C07 and 
whether they describe the same crystal structure with 
various Co contents, is still open. 

In the present work, the basic unit cell, presented in 
Fig. [U with X = 0.50 was employed to account for the 
Y4C03 case. This crystal model has allowed us to use 
the full potential Korringa-Kohn-Rostoker code adapted 
to electronic structure calculations of ordered compounds 
(see. Section nil . 

As far as the superconducting properties are con- 
cerned, it is generally believed, that superconductivity 
in Y4C03 has a singlet-like character and probably medi- 
ated by the electron-phonon interactions The system 
exhibits all the characteristic features of conventional 
superconductors (see Refs. HQ and references therein), 
e.g. perfect diamagnetic Meissner state, the specific heat 
jump at Ts correlated with drop in resistivity and suscep- 
tibility. More puzzling seems the type of its magnetism. 
The system shows e.g. a modified Curie- Weiss behavior, 
with susceptibility maximum — and specific heat jump 
near the Curie temperature,"^^ suppression of magnetism 
with pressure!? and a hysteresis loop in the magnetiza- 
tion measurements^ interpreted as fingerprints of a fer- 
romagnetic state. There are different values of magneti- 
zation reported in literature. For the 4:3 composition one 
can find M = 0.045 ^jlb/^-v.- measured for Y4C03.03 sam- 
ple^ and extrapolated to T = OK or M ~ 0.03 /^^/(f.u. 
Y4C03) at 1.24 K.^ Note, that these values are not exact, 
firstly diamagnetism induced by superconductivity below 
2.5 K does not allow to accurately measure bulk magne- 
tization, and secondly it is difficult to measure magneti- 
zation in weak ferromagnetic systemi^ 

From the NMR measurements^^, it was found that 
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there were three inequivalent Co sites in Y4C03, but 
magnetism was attributed only to the cobalt atoms on 
(2b) site. The estimated magnetic moment was /ico(3) = 
0.23 /is, and important spin fluctuations were observed. 
Actually, there are no other experimental papers on de- 
termination of the Co(3) magnetic moment, thus our 
theoretical result is compared to this value. Most of 
the experimental results concerning magnetism of the 
Y-Co system were interpreted in the framework of the 
spin-fluctuation theory of the weak itinerant ferromag- 
netism,^'^ showing characteristic features as e^. the zero- 
field magnetization M(jf oc T^/^ (Ref. Hj or resis- 
tivity p (X T'^ (Ref. llTl ). However, more exotic mod- 
els without long-range ferromagnetic ordering were also 
considered, i.e. basing on the /iSR measurements^^ the 
magnetic state of Y9C07 was called 'crypto- itinerant fer- 
romagnet', some hybrid magnetic and superconducting 
state was suggested)^ while short range order was dis- 
cussed elsewherei^ 



On the whole, magnetic properties of the Y-Co system 
near 4:3 or 9:7 stoichiometry were found to be similar and 
weakly dependent on Y:Co compositior^ (comparing to 
superconducting behaviors), maybe except for increas- 
ing magnetization with increasing Co content 1^ All these 
data allow to suppose that general conclusions deduced 
from the first principles calculations for Y4C03 should 
also be maintained for Y9C07 (e.g. we may expect sim- 
ilar local magnetic moments). Also, some experimental 
results available only for Y9C07 (e.g. effect of pressure 
on magnetism) , should provide reliable comparison with 
theoretical results obtained for the quasi-ordered model 

of Y4C03. 



To summarize this short review, the commonly ac- 
cepted phenomenological model of the coexistence of su- 
perconductivity and magnetism in Y4C03 / Y9C07 sys- 
tem is as follows: there is an itinerant weak ferromagnetic 
state below Tc with the Co (3) sublattice responsible for 
magnetism, combined with superconductivity below T^, 
being mostly attributed to the Y-Co triangular prisms in- 
side the unit cell. Both physical phenomena visibly com- 
pete (e.g. external pressure suppresses magnetism and 
enhance superconductivity^^"--), but their coexistence is 
possible due to some spatial 'separation' of atom sub- 
lattices responsible for different phenomenai^'i However, 
this separation cannot be treated in the strict sense, since 
the Ginzburg-Landau coherence length was estimated to 
be as large as 300 A (roughly thirty times larger than 
the unit cell size) 1^ One can tentatively imagine this un- 
common state as a 'superconducting sea' with the Co(3) 
'magnetic islands' embedded. The question on the exis- 
tence of the 'magnetic islands' and whether Co(3) mo- 
ments are ordered in Y4C03 has not been yet addressed 
to first principles calculations. Hence, the KKR results 
presented here may help to enlighten this problem. 



II. COMPUTATIONAL DETAILS 

Electronic structure calculations were performed us- 
ing the full potential Korringa-Kohn-Rostoker (FP-KKR) 
technique based on the Green function multiple scatter- 
ing theory?22r— In our implementation of the FP-KKR 
method, the unit cell is divided into the set of gener- 
alized Voronoi polyhedrons arbitrary formed around in- 
equivalent sites, that completely fill the Wigner-Seitz 
cell. The crystal potential has been constructed in 
the framework of the local spin density approximation 
(LSDA) , using Perdew-Wang formulai^ for the exchange- 
correlation part. All results were carefully checked for 
the k-point number convergence, using more than 100 
points in the irreducible part of Brillouin zone for the 
tetrahedron integration method. The results of FP-KKR 
semi-relativistic calculations are presented, with the fully 
relativistic treatment of core electron energy levels. 

In this work, the ordered model of the hexagonal cell 
(Fig. [1]) was used in the electronic structure calculations. 
To take into account the partial filling of the Co(2b) 
site in Y4C03 in KKR method, we employed a 'quasi- 
disordered' structure, i.e the (2b) position was split into 
two non-equivalent sites, (la) and (lb). The (la), i.e. 
(0,0,0) position, was occupied by Co atom (Co(3)) and 
(lb), (0,0, i) site, by a vacancy (Vac), i.e. an empty 
sphere with Z — 0. This structural modification resulted 
in the lowering of the hexagonal cell symmetry from space 
group P63/m (No. 173) to P-3 (No. 147), with reduc- 
tion of symmetry operations from 12 to 6. As above- 
mentioned, this model admits to consider the Y4C03 
composition only. It was verified by additional calcula- 
tions that the KKR results remain unchanged when the 
positions of Co(3) and Vac were exchanged (cobalt on 
(lb) site and empty spheres on (la) site), as one expects 
from equivalence of both sites in the proper space group. 




-0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 

E(Ry) 

FIG. 2: (Color online) Total and site-decomposed densities of 
electronic states in Y4C03 (per unit cell, i.e. three formula 
units), from non-spin-polarized calculations. 



4 



III. RESULTS AND DISCUSSION 
A. Non spin-polarized electronic structure 

The discussion of electronic structure of Y4C03 is 
started from analysis of general features in non spin- 
polarized case. The total and site-decomposed electronic 
densities of states (DOS) are presented in Fig. [2] Due to 
the large number of atoms in the unit cell and in con- 
sequence various interatomic distances, the DOS has a 
complex shape. The Fermi level (Ep) is located in the 
DOS valley, which can tentatively support the chemi- 
cal stability of this system. The valence band region is 
formed from the strongly hybridized atomic states of Co 
(3d, 4s) and Y (4d, 5s), with important s — p and s — d 
charge transfers. The Y-4p states constitute semi-core 
level, located about -1.5 Ry below the Fermi level (not 
shown). Generally, the body of the electronic structure 
is dominated by contributions from the Y(l)-Y(2)-Co(l)- 
Co(2) block, which form the triangular prisms inside the 
unit cell (Fig. [1]). Nominally, the major contribution to 
DOS comes from six Co(l) atoms, due to the highest 
number of valence electrons (54 e) given to the bands. 
The single Co(3) atom, building the separated cobalt 
chains, roughly governs the position of the Fermi level 
due to DOS significantly different from other atoms (a 
high d-like peak below Ep), what enables the Co(3) to 
play important role in ground state properties of Y4C03. 

To have a deeper insight into influence of each atom on 
electronic structure, the site-decomposed DOS are plot- 
ted for inequivalent sites (Fig. |3] and Fig. |4]). The most 
striking DOS feature of Y4C03 is that except for Co(3), 
the Fermi level is systematically placed in the local DOS 
valleys, resulting in relatively low total DOS. The calcu- 
lated small DOSs at Ep per atom (see. Table HIl) would 
rather not suggest a transition to ferromagnetic state. In 
contrast to other atoms, Co(3) exhibits apparently dif- 



TABLE II: Electronic properties of Y4C03. Values of DOS 
n{EF) are given in (Ry~^/spin), magnetic moments in (/is), 
S stands for the Stoner parameter. 



non-spin-polarized calculations 



Atom 


n{EF) 




nd{EF) 


S 


Y(l) 


3.62 




2.50 


0.15 


Y(2) 


4.85 




3.37 


0.19 


Co(l) 


7.71 




5.79 


0.43 


Co(2) 


8.00 




6.86 


0.48 


Co(3) 


10.31 




8.19 


0.61 


Vac 


1.73 




0.47 






N[Ef) 


— 250 Ry ^ per unit cell 




spin-polarized calculations 


Atom 


n-^iEp) nilEp) 


n^diEp) nid{Ep) 




Y(l) 


4.00 


4.22 


2.72 2.91 


-0.024 


Y(2) 


2.90 


4.02 


1.96 2.92 


0.004 


Co(l) 


6.37 


7.18 


4.83 5.18 


-0.007 


Co(2) 


6.45 


7.38 


5.49 6.41 


-0.017 


Co(3) 


4.26 


10.75 


2.37 9.14 


0.551 


Vac 


1.07 


1.15 


0.39 0.37 


0.016 




N^{Ef) = 99 Ry"^ 


N^{Ef) = 120 Ry-i 





ferent DOS, due to a large d-like peak found just below 
the Fermi level (Ep is placed on the decreasing DOS 
slope). This feature yields the highest DOS value (per 
atom) for Co(3) in comparison to other contributions 
(Table HH. Actually, the value n{Ep) ~ 10 Ry-Vspin 
is too small to satisfy the Stoner criterion, since the cal- 
culated Stoner parameter for Co(3) atom S = Idnd{Ep) 
is only about 0.6. Thus, the Stoner analysis based on 
non-spin-polarized DOS and the exchange integral pre- 
dicts a non-magnetic ground state of Y4C03. However, 
the presence of large and narrow DOS peak in the vicin- 




■ Co(l) il 










FIG. 3: Site decomposed DOS per atom from non-spin- 
polarized calculations. 
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FIG. 4: (Color online) Comparison of the Co(3) partial DOS 
from non-spin-polarized (left) and spin-polarized (right) cal- 
culations. Peak under Ef in non-magnetic state is well visi- 
ble. 
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ity of Ep causes that the simple Stoner criterion of ferro- 
magnetism onset should be taken with care and accurate 
spin-polarized calculations may determine the preferred 
ground state. 



B. Ferromagnetism 

Indeed, the spin-polarized calculations do not con- 
firm the non-magnetic state deduced from the Stoner 
criterion. Figure [5] presents the spin-polarized DOS 
shape of Y4C03. At first glance, the differences be- 
tween the non spin-polarized (Fig. [2]) and spin-polarized 
(Fig. [S]) total DOSs are hardly visible, since spin-up and 
spin-down DOS functions are very similar. However, 
the spin-polarized KKR calculations finally resulted in 
stable, ferromagnetic ground state. The partial, site- 
decomposed DOSs (Fig. [5]) show that only the Co(3) 
atom exhibits polarization in DOS shape, as depicted 
in Fig. |4l being confirmed by appearance of local mag- 
netic moment (Table lll| . The Co(3) magnetic moment 
of about /J.C03 = 0.55 fiB decides in favor of magnetism 
in Y4C03, since the remaining atoms possess small mag- 
netic moments (< 0.02 /is). Such a small DOS polar- 
izations seen on other atoms should be rather considered 
as a response to magnetic field of Co (3) atom, than an 
intrinsic local magnetic moments. The total magnetiza- 
tion Mtot — 0.38 ^s/unit cell yields M ~ 0.13 fiB per 
Y4C03 formula unit. Mtot is lower than ^co3 due to 
the overall 'diamagnetic' response of the Y-Co triangular 
prisms, generating about -0.1 fj,B- The comparison with 
the experimental values is discussed in Section IIII Dl 

The specific character of the Y4C03 electronic struc- 
ture can be observed in the electron dispersion curves, 
shown in the Fig. [HI In the T — M — K — T triangle 
no band crosses the Fermi level for both spin directions, 
and this behavior holds for any direction in the — 
plane, as was checked by extensive calculations along ran- 
domly chosen directions. Similar result was earlier ob- 
served from non spin-polarized calculations The fer- 
romagnetic state leads mainly to the slight shift of dis- 
persion curves and some changes in details. Since, the 
Y4C03 system is predicted to possess an energy band 
gap in the fcz = plane, it would be interesting to verify 
it (likely detectable on the ARUPS spectra, but due to 
lack of single crystals there are no such results available). 
The vanishing Fermi surface in the kz = plane is pre- 
sumably responsible for the DOS valley formed around 
Ep (Fig.[5l). In the T — A direction (along the k^ axis), al- 
most linear and strongly dispersive bands cross Ep, while 
in the A — L — H — A triangle, which is just a shift of 
the r — M — K — T triangle to the Brillouin zone top, 
bands crossing Ep become flat, which mostly gives rise 
to metallic- like properties of Y4C03. 

On the whole, the important result from presented 
spin-polarized KKR calculations, concerns the fact that 
there is no need to search for unconventional mechanism 
to explain magnetism in Y4C03 system. The energeti- 
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M,„, = 0.38 /(B/unit cell 
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FIG. 5: (Color online) Total and site-decomposed densities of 
electronic states in Y4C03 from spin-polarized calculations. 



cally unfavorable DOS peak just below Ep (Fig. \^ on 
the single Co(3) atom, seems to be the reason of turn- 
ing the system into the ferromagnetic state. In magnetic 
state, the spin- up peak is not seen near Ep, while in 
the spin-down DOS the peak is smaller and is expelled 
above Ep . In Ref . [Ill, devoted to the electronic structure 
and magnetism of Y9C07, authors suggested two mod- 
els of magnetism of this system, called A and B. Briefly, 
in the model A magnetism was attributed to the excess 
Co atoms on Y(l) and Y(2) sites, while in the model B 
the magnetism came from itinerant electrons of Co on 
(2b) sites. Authorsii did not definitely conclude, which 
model was valid, but the model A was said to be the most 
probable, and the 'perfect Y9C07 crystal' was expected 
to be paramagnetic. From our results we see, that itin- 
erant ferromagnetism is an intrinsic property of Y4C03 
system. Hence, additional defects as the excess Co atoms 
on Y sites are not necessary to give rise ferromagnetism, 
even if instead of long range order in real material short 
range order appears, as suggested e.g. in Ref. [2^. Simi- 
lar behavior is expected to be valid for YgCoy, since as 
we mentioned in Introduction, experimentally the mag- 
netic properties of the Y-Co system do not change much 
between Y4C03 and Y9C07.— 

To look closer on the magnetism of the Co(3) atom we 
can calculate the spin dependence of d orbital occupation 
(actually, the diagonal elements of occupation matrix cal- 
culated using imaginary part of Green's function, defined 
as (/imil — ilmG'|^2W2)). The Co(3) (2b) site, being 
surrounded by Y(l) triangles in z = ±|; planes, has a 
trigonal, C^i symmetry. Thus, there are two degenerated 
Eg doublets, {dxy,d^2_y2} and {dyz,dxz\ with occupa- 
tion for both doublets approximately equal to O.Sle-^ and 
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0.76e4_. The dz2 orbital {Ag representation, Z = 2, m = 0), 
aligned in the z direction, is the highest-energy and least 
occupied orbital. It also has the highest level of polariza- 
tion, with 0.79e^ and 0.49e^. Hence, about 60% of mag- 
netization conies from d^^ orbital and microscopically it 
is the main source of magnetism in this compound. The 
directional character of this orbital favors the simple fer- 
romagnetic ordering of magnetic moments along the z- 
axis. 

The fact that the magnetic moments are 'localized' 
only on Co(3) atoms can be well visualized on the spin 
magnetization contour maps, presented in Figs. [7]and|Sl 
The spin magnetization is calculated as a difference be- 
tween the spin-up and spin-down electron densities. The 
magnetic moments resides on the Co(3) corner atoms, 
while the rest of the presented xz plane has magnetiza- 
tion close to zero (Fig. [S]). The shape of magnetization 
around the Co(3) atom in the xz plane (Fig. [7^) confirms 
that the 3d^2 orbital plays the major role in constitution 
of magnetic state. The xy cross-sections well refiect the 
local symmetry of the magnetization, which should be 
trigonal (Csi). Interestingly, in the z = plane we may 
have an impression of cylindrical symmetry (Fig.[7]3), due 
to the equally distanced Y(l) triangles (z = ij), each 
rotated by 60°. But, when we move off the z = plane, 
the triangle symmetry occurs (Fig. [Tj;). 

Magnetism in Y4C03 appears to be unusual. There 
is only one 'magnetic atom' in the complex unit cell, 
and twenty other atoms form diamagnetically polarized 
background. The reason why other Co atoms are non- 
magnetic is understood from non spin-polarized partial 
DOS, where only Co(3) atom exhibits unstable DOS 
function due to the pronounced peak near Ep, while for 
Co(l) and Co(2) Ep is located in the DOS minimum. In- 
ducing magnetic moments on the Co(l) and Co(2) atoms 
would make the system energetically unfavorable. 

Another peculiarity of this crystal structure is the ar- 
rangement of Co(3) atoms in chains, lying in channels 
formed by the octahedras of Y(l) atoms (see. Fig [T|). 
Thus, the Co(3) atoms form a quasi-one-dimensional 
magnetic structure, a feature which was announced for 
this system2i!^ but not explored in the past. This low- 



dimensionality rises the question if the long-range ferro- 
magnetic order can occur here, as the Mermin- Wagner 
theorem'^^ predicts the lack of long-range magnetic or- 
der within a Heisenberg model with short-range interac- 
tions in one dimension. This of course does not com- 
pletely rule out the presence of magnetic ordering for the 
three dimensional structures with quasi-one-dimensional 
chains, but the properties of such systems are affected by 
the low-dimensional character. Recently, the quasi-one- 
dimensional ferromagnetism was investigated in some 
cobalt-based systems, like monoatomic Co chains on Ft 
substratc'^^ or Co oxides, like BaCoOs or Ca3Co206 
(see, e.g. Refs. I38l - |40[) . It is interesting to underline 
some structural similarities between Y4C03 and afore- 
mentioned Co oxides. The Co(3) chains surrounded by 
Y(l) triangles, forming octahedras in 2;-direction, are 
similar to CoOg octahedras in oxides family. In both 
cases, the chains are running along c axis of the hexago- 
nal cell. The geometrical difference is that in Y4C03 the 
intra- and inter-chain Co-Co distances are about twice as 
large as in the oxides. The large interchain distance in 
Y4C03 (11.5 A) rather excludes the importance of inter- 
chain magnetic interaction, which is of great importance 
in oxides, and leads to two-dimensional antiferromagnetic 
effects.-SS. This quasi-one-dimensional character of mag- 
netism in Y4C03 makes this system even more interesting 
and opens new areas for future investigation. 



C. Ferromagnetism vs. superconductivity 

The fact that most of the unit cell volume of Y4C03 
has negligible magnetization, allows to draw some qual- 
itative conclusions on the possibility of coexistence of 
ferromagnetism and superconductivity. For example, in 
the xz cross-section of the unit cell, presented in the 
Fig. m about 90% of the area has magnetization lower 
than lO^'^ ^B/a%. For the whole unit cell, only 1.3 % of 
the volume has magnetization higher than 10^^ ^b/clb- 
Roughly speaking, if the electron density inside the unit 
cell has negligible polarization, in principle it can not dis- 
turb forming singlet Cooper pairs. Certainly, the Cooper 
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FIG. 8: (Color online) Spin magnetization distribution in the 
y = plane (face of the unit cell), showing that the charge 
density in the large areas of the unit cell is almost not polar- 
ized. 



pairing acts in momentum space, so we should gener- 
ally verify whether there are (k-|-, — kj,) electrons near 
the Fermi level. Since the (k, — k) degeneracy is en- 
sured by the centrosymmetry of the unit cell, we should 
look for the (k-t-,k_i.) states. From the band structure 
plot e.g. in the T — A direction we observe three bands 
crossing Ep- Two of these bands are almost the same 
for both spin-directions and they are crossing Ep in 
points: k^ = (0,0,0.2308)^, k^ = (0,0,0.2307)^ and 
kt = (0,0,0.3422)^, k_^ = (0,0,0.3410)^. These bands 
are mostly of the Co(l) and Y(2) character. So, one 
concludes that in principle there are electrons in the Y- 
Co trigonal prisms, which can form Cooper pairs and 
lead to singlet superconductivity. But more quantitative 
analysis of electron-phonon coupling is required to assess 
whether the BCS type superconductivity can appear in 
Y4C03. 

Nevertheless, our results and conclusions may support 
the model of coexistence of weak ferromagnetism and su- 
perconductivity in Y4C03, with the ferromagnetism car- 
ried by chains of Co(3) atoms, screened by superconduct- 
ing trigonal prisms of Y(l)-Co(l)-Y(2)-Co(2). Interest- 
ingly, this picture is qualitatively similar to the vortex 
lattice in type-II superconductors, since in both cases we 
have superconducting sample penetrated by the magnetic 



field lines, forming hexagonal lattice. 



D. Spin fluctuations 

The calculated value of Y4C03 magnetization 
(0.13 //^/(f-u.)), if compared to the measured 
M = 0.045 /^^/(f.u.) is about 2.5 times overesti- 
mated. Similarly, comparing the calculated Co(3) 
magnetic moment /icos = 0.55 fiB with the NMR 
estimations^ HCo^ = 0.23 fiB we get 2.5 times over- 
estimation. This shows, that Y4C03 can be a rare 
example of weak ferromagnetic system, where LDA 
tends to overestimate the tendency to magnetism, and 
suggests that it may be near the ferromagnetic quantum 
critical point.— The well-known similar examples are: 
ZrZna Scgln and Ni3Al/Ni3Ga (see, e.g. Refs. 
where due to strong spin fluctuations appearing in real 
samples, measured magnetic moments are much smaller 
than LDA (or GGA) values. Thus, our KKR-LDA 
results strongly support the classification of Y4C03 as a 
weak itinerant ferromagnet with spin fluctuation effects, 
as already suggested, basing on the analysis of experi- 
mental results. 1^ In Y4C03, overestimation of magnetic 
moment is a little smaller, comparing to aforementioned 
cases, since in hexagonal Sc3ln LDA magnetic moment 
fiSc — 0.35 fiB is about 7 times higher than experimental 
one^^i^ while the factor about 3 is established in ZrZn2 
(Ref. Hi) and Ni3Al (Ref. IH). This suggests that the 
spin fluctuations in Y4C03 should be comparably weaker. 
The spin fluctuations strength parameter Asf, which can 
be treated as the analog of the electron-phonon coupling 
constant Aph, can be extracted from the electronic 
specific heat coefficient 7. The experimental value lies 
in the range of 7 ~ 3.1 - 3.4 mJ/(mol at. K^). 2SJ^ 
Using the KKR-LDA value of total DOS at Fermi 
level (spin-polarized case), N{Ef) — 229 Ry~^ we 
get 7LDA — 1-89 (mJ/mol at. K^). Assuming that 
the band value is renormalized by the electron-phonon 
as well as the electron-paramagnon interaction (spin 
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FIG. 9: Evolution of the Co(3) magnetic moment witii pres- 
sure. 

fluctuations), the formula 7 — 7lda(1 + Aph + Asf), 
gives Aph + Asf = 0.6 — 0.8. This allows to safely put 
the upper limit Aph + Asf < 1. It is possible to have 
independent estimation of Aph and Asf if one admits 
that superconductivity is driven by phonons and the 
McMillan formula^* can be applied here. We use the 
experimental value of and the renormalized McMillan 
formula to take into account the spin fluctuations)^ 
(see, also Ref. 50: and references therein) 

J. _Qd [ 1-04(1 + Aeff) ] 

" 1.45 ""^1 AefF-MefF*(l + 0.62AefF)J ' ^' 

with the renormalization: 

Aeff = Aph/(1 + Asf) 
Aieff" = (/X*+Asf)/(1+Asf). (2) 

The experimental Debye temperature is Qd — 215 K.^* 
Noteworthy, the strong renormalization of Coulomb 
pseudopotential parameter fi* restricts the possible range 
of Asf, e.g. we can get ~ 2.5 K for Asf — 0.1, Aph = 0.7 
and = 0.08 (yielding /ieff* = 0.164 and Aph + Asf = 0.8) 
or for Asf = 0.15, Aph = 0.85 (^eff* — 0.20 and Aph + Asf = 
1.0). Thus, we can estimate the electron-paramagnon 
interaction parameter to be of the order of Asf 0.1, 
but rather not higher than 0.2. This qualitative discus- 
sion shows that the spin fluctuations are likely present 
in Y4C03, but they are not as strong as e.g. in Scain, 
where Asf > 1 can be deduced from the specific heat mea- 
surements in the magnetic fieldi ^^i^^ Noteworthy, similar 
measurements, i.e. the influence of high magnetic field 
on the electronic specific heat, would be very helpful to 
study the presence and strength of spin fiuctuations in 
Y4C03. 

E. Effect of pressure 

The effect of pressure on the magnetic properties of 
Y4C03 was also analyzed. We have performed calcula- 



tions under 'quasi-hydrostatic' conditions. The a and 
c parameters were equally contracted in the range 0%- 
2% and all self-consistent calculations were started from 
the beginning. The bulk modulus, estimated under these 
conditions, is about B — 75 ±5 GPa, which is higher than 
the bulk modulus of yttrium (^ 40 GPa^) and much 
smaller than the corresponding value for cobalt 200 
GPa—). In calculations for the smaller unit cell volume, 
the magnetism is suppressed and the magnetic moment 
decreases rapidly, with the slope ^ ~ 0.08 ^s/GPa, 
as seen in Fig. [HI The reason for such behavior can be 
apparently seen from the evolution of the partial Co(3) 
DOS under pressure fPig. fTO]) . When the unit cell volume 
decrease, the spin-up DOS peak under Ep tends to con- 
stitute, while the spin-down peak, expelled above Ep in 
magnetic state, tends to move below Ep. Thus, the pres- 
sure lowers the DOS polarization and the magnetic mo- 
ment, which additionally confirms that the Co(3)-DOS 
peak is responsible for appearance of magnetic moment 
on cobalt atoms on (2b) sites. 

The suppression of magnetism under pressure from 
KKR calculations remains in agreement with the experi- 
mental trends observed in Y9C07 (e.g. lowering of Curie 
temperature and magnetizatioufiSi^ for Y4C03 there are 
no such results available in the literature) . Extrapolation 
of the curve from Fig. [5] to zero magnetic moment gives 
the critical pressure, where magnetism completely disap- 
pears, as Pc — 7 GPa. This value can be compared with 
very recent measurements of Y9C07 under pressure^, 
where the critical pressure was estimated to be Pc — 3 
GPa. The overestimation of the critical pressure is in line 
with the spin fluctuations and weak ferromagnetism in 
this system, as revealed from our calculations. Interest- 
ingly, if the zero-pressure value of the magnetic moment 
the NMR experimental estimation /j,(Co3) ^ 0.23 fiB is 
accounted for, the theoretical slope ^ ~ 0.08 ^s/GPa 
also predicts decrease of ^ to zero for p ~ 3 GPa. 



IV. SUMMARY AND CONCLUSIONS 

The results of the FP-KKR electronic structure cal- 
culations for Y4C03 system were presented. The ferro- 
magnetic state obtained from spin-polarized computa- 
tions can be attributed to the single Co atom located 
on the (2b) site, being the only magnetic atom among 
21 ones in the unit cell, and forming the quasi-one- 
dimensional magnetic chains. The LDA values of the 
magnetization (M ~ 0.13/i_B/f.u.), Co(3) magnetic mo- 
ment (fi ~ 0.55/xb) and critical pressure {pc ~ 7 GPa) are 
overestimated comparing with experiments, which can be 
tentativelly explained in terms of weak ferromagnetism 
with moderate spin fluctuations (Asf ~ 0.1). The cal- 
culated spin magnetization distribution as well as other 
band structure parameters, tend to the conclusion that 
the conventional, singlet-like superconductivity may co- 
exist with ferromagnetism in Y4C03, due to relatively 
weak magnetic moments arranged along thin chains (the 
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unit cell edges) on the one hand, and the presence of 
non-polarized electrons at the Fermi level (filling most of 
the unit cell) on the other hand. On the whole, the FP- 
KKR results confirmed that the quasi-one-dimensional 
magnetism is an intrinsic property of Y4C03 and its co- 
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